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Abstract. We study the temperature dependence of bottomonium for 
temperatures in the range 0.4rc < T < 2.1Tc, using non-relativistic 
dynamics for the bottom quark and full relativistic lattice QCD simulations 
for Nj = 2 light flavors. We consider the behaviour of the correlators in 
Euclidean space, we analyze the associated spectral functions and we study 
the dependence on the momentum. Our results are amenable to a successful 
comparison with effective field theories. They help build a coherent picture 
of the behaviour of bottomonium in the plasma, consistent which the current 
LHC results. 



1. Introduction 

In this note we review the results on the spectrum of bottomonium in the quark- 
gluon plasma obtained by our FASTSUM collaboration [T]-[3] using a hybrid 
approach which combines exact relativistic dynamics for the light quarks with 
a non-relativistic approach for the bottom quark. 

Recently, the first detailed results on the sequential Upsilon suppression in 
PbPb collision at LHC energies have appeared 0] (see also several talks at 
this meeting [5]). These results call for a quantitative understanding of the 
behaviour of bottomonium at high temperature. Due to their considerable 
experimental and theoretical interest these studies have been carried out in the 
past with a variety of techniques. We refer to recent excellent reviews [SlEj, as 
well as to our published papers, for background material and a complete set 
of references to early studies; here we limit ourselves to an exposition of our 
approach and results. 

In the next Section we outline our approach to the problem. Next, we present 
the information we have obtained from the correlators alone. We continue by 
showing the results on the spectral functions. We finally discuss momentum 
dependence and we close with a short summary. 

2. Five steps to Bottomonium 

Our strategy is to treat as accurately as possible the light quarks, and use non- 
relativistic QCD (NRQCD) for the b quarks. In this way we hope to keep the 
systematic effects under control. Moreover, the accuracy of the results can be 
improved following the standard lattice techniques. In short outline, this is our 
'five steps approach to bottomonium' strategy: 

i) Use full relativistic dynamics for up and down quarks (the inclusion of 
the strange quark is in progress). Gauge configurations with two degenerate 
dynamical light Wilson-type quark flavors are produced on highly anisotropic 
lattices, in a range of temperatures comprised between OATc and 2.1Tc. Details 
of the lattice action and parameters can be found in Refs. [51[S]. 

ii) We use NRQCD for the bottom quarks. Since NRQCD relies on the scale 
separation M ^ T and we study temperatures up to 2Tc ~ 400 MeV, its 
application in our temperature range is fully justified. The light quarks, instead, 
feel the full relativistic dynamics. We computed NRQCD propagators on our 
configurations using a mean-field improved action with tree-level coefficients, 
which includes terms up to and including 0{v'^), where v is the typical velocity 
of a bottom quark in bottomonium. 

iii) Analyze bottomonium correlators in Euclidean space. In standard 
conditions we expect an exponential decay. An analytic model for the Euclidean 
thermal behaviour is not known — we contented ourselves to study what to 
expect when quarks are no longer bound. Consider free quarks in continuum 
NRQCD with energy Ep = p^/2M. The correlators for the S and P waves are 




Figure 1. Effective mass plots for the T (left) and Xbi (right) using point 

sources for various temperatures: note the different temperature dependence 
(Ref. m). 

then of the form [12] 

Gsir) ~ / exp(-2i?pr) ~ r"^/^, (i) 

Gp(t) ~ / p' exp(-2i?pr) ~ r^^/^, (2) 

i.e., they decay as a power for large euclidean time. Of course, interactions 
and finite lattice spacing and volume effects are expected to modify this in the 
realistic case. Moreover, again, this is a limiting high-temperature behaviour 
which should not be expected to hold even if the quarks unbind in the plasma. 

iv) Obtain spectral functions. In NRQCD, the spectral relation reduces to 

exp(-a;'r)/)(w') (NRQCD), (3) 

even at nonzero temperature (note w = 2M + co'). As a result, all problems 
associated with thermal boundary conditions are absent. We will extract 
NRQCD spectral functions from the euclidean correlators using the Maximum 
Entropy Method (MEM) [T7] . 

v) Extract physics from the results, discuss phenomenology, validate effective 
models. We underscore that any meaningful comparison with experimental 
results and continuum physics requires control over systematic and lattice 
artifacts. 



3. T and Xb correlators in the plasma 

We discuss here the results one can obtain from the analysis of correlators in 
Euclidean time, focusing on the T and Xb states [1]. 

In Fig. 1, standard effective masses, defined by mefj(T) = — log[G(T)/G(r — 
ttr)], are shown for both the T and the Xbi propagators at various temperatures. 
Single exponential decay should yield a r-independent plateau. In both cases 
we find that at the lowest temperature, T = 0.42Tc, exponential behaviour is 
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Figure 2. Effective exponents 7efr('7") for the T (left) and Xbi (right), as a 
function of Euchdean time for various temperatures. The dotted hne indicates 
the non-interacting result in the continuum (Ref. [T]). 



visible provided one goes to late euclidean times. At higher temperature we can 
see first indications that the T is not sensitive to the quark-gluon plasma up 
to T ~ 2Tc, while the Xb may melt at much lower temperatures, as exponential 
behaviour is no longer visible. 

To visualize the approach to quasi-free behaviour we construct effective 
power plots jT|, using the definition 7eff(r) = —7"^^ = — r '^'-'^'''"^^ J^j:^ , 
where the prime denotes the (discretized) derivative. For a power decay, 
G{t) ~ T~'^ , this yields a constant result, 7efr(T) = 7. The results are shown 
in Fig. 2. We confirm again that the T displays a very mild temperature 
dependence, while for the Xbi the effective power tends to flatten out. As 
discussed, in the continuum limit and in the free case one should observe a flat 
behaviour. The observed flattening might thus suggest an approach to a free 
field behaviour, however several caveats apply. Also shown are the effective 
exponents in the continuum non-interacting limit. In the case of the Xbi^ we 
observe that the effective exponent tends towards the non-interacting result at 
the highest temperature we consider. 

To extend and validate these results, we proceed to the calculation of the 
spectral functions. 



4. Spectral functions 

Our main results for the spectral functions [2 can be seen in Fig. [3l which 
shows that as the temperature is increased the ground state peaks of the T and 
rjh remain visible. One caveat applies also here: the apparent width at zero 
temperature is most likely due to a lattice/MEM artifact, and this calls for an 
analysis of the discretization effects on the spectral function, which is one of our 
ongoing projects. The peaks associated to the excited states become suppressed 
at higher temperature and are no longer discernible at T/Tc ~ 1.68. The 
temperature dependence of the position and width of the ground state peaks 
is compared to analytical predictions obtained within the EFT formalism |18j . 
We note that the survival of the T(IS') state and suppression of excited states 
is consistent with the recent experimental results 
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Figure 3. Spectral functions p{uj), normalised with the heavy quark mass, 
in the vector (T) channel (upper panel) and in the pseudoscalar (77^) channel 
(lower panel) for all temperature available. The subpanels are ordered from 
cold (top left) to hot (bottom right). Every subpanel contains two adjacent 
temperatures to facilitate the comparison (Ref. [2]). 

From these spectral functions we can determine the mass (from the peak 
position) and (an upper bound on) the width of the ground state at each 
temperature. In Fig. S] we show the temperature dependence of the mass shift 
A£^, normalised by the heavy quark mass and the temperature dependence of 
the width, normalised by the temperature. 

We now contrast our results with analytic predictions derived assuming a 
weakly coupled plasma. According to Ref. |18j , the thermal contribution to the 
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Figure 4. Position of the ground state peak Ai?, normalised by the heavy 
quark mass (left), and the upper limit on the width of the ground state peak, 
normalised by the temperature (right), as a function of T/Tc in the vector (T) 
channel. Similar results have been obtained for the pseudoscalar (?];,) channel 
(Ref. m- 



width is given, at leading order in the weak coupling and large mass expansion, 

by 

r 1156 



^ ^^'-al^U.27al (4) 

i.e., the width increases linearly with the temperature. If we take as an estimate 
from our results that T/T ~ 1, we find that this corresponds to ~ 0.4, which 
is a reasonable result. (It would be of interest to compute as directly on our 
configurations.) In the same spirit the thermal mass shift is given in Ref. [18] 
by 

(5-Bthermal = Cts — 5.93as . (5) 

Lnermai 9 M M 

In these simulations we have Tc ~ 220 MeV, M ~ 5 GeV. Taking these values 
together with as ~ 0.4 as determined above, Eq. ([5]) becomes 



M \M J \TJ \Tcy 

In order to contrast our results with this analytical prediction, we have 
compared the temperature dependence of the peak positions to the simple 
expression 

AE . 



M + (f;j . P) 

where c is a free parameter. This is shown by the dashed line in Fig. U] (left 
panel). The numerical results and the analytic ones are not inconsistent, within 
the large errors, which encourages us to pursue further these studies. 



5. Momentum Dependence 

We extend now our study of bottomonium spectral functions in the quark-gluon 
plasma to nonzero momentum [3]. Note that effective theories analysis predict 
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Table 1. Nonzero momenta used in this study. Also indicated are the 
corresponding velocities v = \p\/Ms of the ground states in the vector (T) 
and pseudoscalar (ry^) channels, using the ground state masses determined 
previously P, Mx = 9.460 GeV and M^^, = 9.438 GeV. 



significant momentum effects at large momenta, and that current CMS results 
have been obtained at large momenta. So there is both phenomenological and 
experimental motivation for these studies. 

The momenta and velocities that are accessible on the lattice are constrained 
by the discretization and the spatial lattice spacing. The lattice dispersion 
relation reads 

2 2 ■ 2Pi 27rni Ns Ns 

a.P =42^sm -, pi = ——, - <ni< —. (8) 

i=i * 

To avoid lattice artifacts, only momenta with rii < Ns/4: are used: we 
consider the combinations (and permutations thereof) given in table [TJ The 
largest momentum, using n = (2,2,0), is |p| ~ 1.73 GeV, corresponding to 
V = \p\/Ms — 0.2. Therefore, the range of velocities we consider is non- 
relativistic. 
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Figure 5. Position of the ground state peak M(p)/M(0) (left) and the upper 
limit on the width of the ground state peak, normalized by the temperature, 
r/T (right), as a function of the velocity squared (f^) in the vector (T) channel 
Analogous results have been obtained for the the pseudoscalar (rn,) channel. 
The dotted line in the left figure represents M(p)/M(0) = 1 + 1^2 ^j^gf_ jg]). 

As in the zero momentum case, we extract masses and widths from our 
correlators using MEM, and we display them in Fig. 5. We observe that the 



peak position increases linearly with f^, as expected. Assuming the lowest- 
order, non-relativistic expression Af (p) = M(0) + p^/2M(0), one finds 



M(0) "'^+2M2(0) "^+2'' ' 

which is indicated with the dotted lines in the left figures. 

The dependence on the velocity can be compared with EFT predictions. In 
Ref. |19j a study of the velocity dependence was carried out in the context 
of QED, working in the rest frame of the bound state (i.e. the heat bath is 
moving). In order to compare with our setup, we consider the case in which the 
temperature is low enough for bound states to be present and that the velocities 
are non-relativistic. In that case, one finds jl9j . in the rest frame of the bound 
state and at leading order in the EFT expansion. 



^,,Jl±l]. (10) 
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where Fq is the width at rest. Interpreting the width as an inverse lifetime, one 
can express this result in the rest frame of the heat bath by dividing with the 
Lorentz factor 7 = — v'^. An expansion for non-relativistic velocities then 

yields 

^ = l-^ + 0(-^). (11) 
i o 

If we apply this result to our study of bottomonium, we find that the effect 
of the nonzero velocity shows up as a correction at the percent level (recall 
that < 0.04), which is beyond our level of precision but consistent with the 
observed v independence within errors. Similarly, additional thermal effects in 
the dispersion relation are currently beyond our level of precision. In summary, 
the observation in our low- momentum range are consistent with Ref. [19], and 
in order to observe the predicted non-trivial momentum dependence we need 
to explore larger momenta. 



6. Summary 

We have presented our results on bottomonium in the quark-gluon plasma, for 
temperatures up to 2.1Tc, at the threshold of the region currently explored 
by LHC heavy- ions experiments. Our analysis uses full relativistic gauge field 
configurations, and a non-relativistic approach for the bottom quarks. Our 
goal is to reach a quantitative, systematic control on the results, which is called 
for by the new level of experimental accuracy. Related to this, we aim at 
understanding in detail the range of validity of the available effective models. 
This is a task of foremost importance, since an accurate analytic modeling is 
needed as an input to real time evolution. 

Our results so far were based on a set of gauge field configurations obtained 
with two flavors of dynamical quarks. We have used anisotropic lattices in 
order to maximize the number of points in Euclidean time — a much needed 
feature for a successful MEM analysis. In this way even on our hottest lattices 



(T = 2.1Tc) we have been able to analyze correlators and proceed to the 
calculation of the spectral functions. 

On these lattice ensembles we have observed that the temperature affects 
only very mildly the T and ijh correlators, suggesting the survival of the 
fundamental states in these channels up to T ~ 2Tc. A subsequent calculation 
of the spectral functions carried out with a MEM analysis has confirmed the 
survival of the fundamental states and revealed the melting of the excited states 
in the same channels above lATc. The correlators for the Xb and indicate a 
quick melting of these states above Tc- 

On the methodological side, we have argued that the peculiar nature of 
non-relativistic dynamics greatly simplifies the MEM analysis and enhances its 
reliability. We have invested a considerable effort in monitoring the systematic 
effects associated with the MEM analysis, which is described in our papers, and 
we feel that we have them well under control. The special form of the NRQCD 
spectral relations shows that the NRQCD spectral function is an inverse Laplace 
transform of the correlator. In a nutshell, this makes the MEM analysis easier 
than in the full relativistic case, and perhaps opens the way to a direct, more 
simple evaluation of the spectral functions themselves. 

Our activities in the near future are evolving along different lines. Firstly, 
new results which will be available in the near future include a MEM analysis 
for the Xb and /if, channels. 

Secondly, we wish to approach more closely the continuum limit, with a 
physical matter content. To this end, simulations with a reduced lattice spacing, 
and the inclusion of the strange quarks, are already in progress. On the 
NRQCD side, on the new Nf = 2 + 1 lattices the quark mass is tuned precisely 
using the kinetic mass and the lattice NRQCD dispersion relation. Moreover 
we are considering using the 1-loop determinations of the a ■ B term in the 
NRQCD coefficients which has recently produced an impressive improvement 
in the determination of the (zero temperature) hyperfine splitting [20] . All this 
will enhance the control over various lattice systematic errors. In a longer term 
perspective, relativistic beauty might also become feasible on these finer lattices, 
mirroring the approach that is already being followed for charmonium [8tl21j. 
This will provide a definitive check of the NRQCD approach. 

Next, as for the momentum dependence, we have mentioned that both 
experimental and phenomenological activities call for calculations at larger 
momenta: a sizable dependence is expected at large momenta and large 
momenta are those studied by CMS. Large momentum studies will be possible 
on our finer lattices and this is being planned. 

As a final remark, leading to further developments, an important aspect 
of our study is the comparison with effective theories. Within our limited 
precision, and limited explored range, the comparison has been satisfactory. 
This motivates us to enhance this analysis in the future. To carry out this 
program, besides the improvement already described, we will have to explore 
a wider set of parameters. One economic way of doing this is to vary the 
heavy quark mass, studying artificial mesons: this has been very helpful at zero 
temperature in cross checking lattice results and effective model predictions — 
for instance the results for the light meson spectrum are being compared with 



chiral perturbation theory and the heavy meson spectrum with heavy quark 
effective theory. In the same spirit we would like to know the dependence 
of the bottomonium spectrum on the heavy quark mass, and temperature, 
and compare it with various effective models in their respective ranges of 
applicability. Remember in fact that different effective models hold true in 
correspondence with different relative values of the relevant scales such as 
M,Mas,T,Ma1,mD, so upon varying mass and temperatures one should 
observe a crossover between different behaviours, thus building confidence 
in using the corresponding models. This program has started and its first, 
preliminary results have been reported at Lattice 2012 [22j . 
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